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Abstract. The thermal emission detected from the millisecond pulsar J0437-4715 is not explained 
by standard cooling models of neutron stars without a heating mechanism. We investigated three 
heating mechanisms controlled by the rotational braking of the pulsar: breaking of the solid crust, 
superfiuid vortex creep, and non-equilibrium reactions ("rotochemical heating"). We find that the 
crust cracking mechanism does not produce detectable heating. Given the dependence of the heating 
mechanisms on spin-down parameters, which leads to different temperatures for different pulsars, 
we study the thermal evolution for two types of pulsars: young, slowly rotating "classical" pulsars 
and old, fast rotating millisecond pulsars (MSPs). We find that the rotochemical heating and vortex 
creep mechanism can be important both for classical pulsars and MSPs. 
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INTRODUCTION 

For all standard cooling models [13], neutron stars (NSs) cool down to surface tempera- 
tures T s < 10 4 K within less than 10 7 yr. However, the observation of ultraviolet thermal 
emission from millisecond pulsar J0437-4715 [8], whose spin-down age is T s d ~ 7 x 10 9 
yr [5], suggests a surface temperature ~ 10 5 K for this pulsar. Thus, a heating mechanism 
has to be added to the thermal evolution models in order to obtain agreement between 
theory and observation. 

The goal of this work is to provide a comparative analysis of the thermal evolution 
including three internal heating mechanisms that are driven by the rotational evolution: 
crust cracking [2, 3], vortex creep [1, 12, 9], and rotochemical heating [11]. The first 
operates in the nuclear lattice that composes the crust of the star; the second, in the 
inner crust, where superfiuid neutrons coexist with the nuclear lattice; and the last one, 
mainly in the core of the star. Due to the dependence of heating mechanisms on spin- 
down parameters, we study separately the thermal evolution for two types of pulsars: 
"classical" pulsars, characterized by strong magnetic fields and relatively slow rotation, 
and MSPs with low surface magnetic fields and fast rotation. 



INTERNAL HEATING MECHANISMS 



Typically, NSs are born with temperatures ~ 10 K. The subsequent cooling is caused 
by neutrino emission from the interior and photon emission from the surface of the star. 
The evolution of the internal temperature is given by the thermal balance equation 

t = ^(L H -L Y -L v ), (1) 

where C is the total heat capacity of the star, L y is the photon luminosity, L v is the 
neutrino luminosity, and L# is the power generated by internal heating mechanisms, of 
which we consider the following three: 

Crust Cracking: The rotational braking causes a gradual change in the shape of the star 
from an ellipsoidal to a spherical form, which increases the stress in the crust. When 
it reaches a critical deformation, it breaks, and part of the accumulated strain energy is 
released. The energy loss rate by this mechanism, based on the formalism of Baym & 
Pines [2], is L cc = 5bI0 c ClCl, where b is the rigidity parameter of the crust, / is the mo- 
ment of inertia of the star, C is the maximum strain angle, and Q. is the angular velocity 
of the NS. With b ~ 1(T 7 [4] and 6 C ~ 1CT 1 [7], we find L cc ~ 10 26 erg s" 1 for typical 
MSPs, too low to produce detectable heating. 

Vortex Creep: When a superfluid is forced to rotate, this generates vortex lines whose 
distribution and dynamics are determined by the rotation rate. In a rotating NS contain- 
ing a neutron superfluid, the vortices interact with the nuclei of the inner crust and are 
dragged through the nuclear lattice due to the spin-down of the star. The pinning and un- 
pinning of the vortex lines with respect to the nuclei of the crystal lattice releases energy 
that heats the star. The energy dissipation rate by this mechanism is given by L vc = JQ. 
[1], where the parameter J contains all the (highly uncertain) information about the 
vortex-nuclei interaction. 

Rotochemical Heating: In chemical equilibrium, in a NS composed of neutrons (n), 
protons (p), and leptons (/: electrons and muons), the chemical potentials /i, satisfy 
r\ np i = \x n — fi p — fa =0. However, if the rotation of the star is slowing down, the cen- 
trifugal force is reduced, the central density of the star increases, and the chemical po- 
tentials are imbalanced, r\ np i ^ 0. As the equilibrium composition is altered, the NS 
will tend to relax to the new chemical equilibrium (via beta and inverse beta decays), 
releasing energy as neutrinos, antineutrinos, and thermal photons. The evolution of the 
chemical imbalance is of the form i] np i = —A{r\ np i,T) —R np iQ.Cl [6], where the func- 
tion A quantifies the effect of reactions toward restoring chemical equilibrium, and the 
constant R np \ quantifies the departure from equilibrium due to the change in the angu- 
lar velocity Q. of the star. Reisenegger [11] found that, if the angular velocity Q. varies 
slowly over the times required to cool the star and achieve chemical equilibrium, the star 
reaches a quasi-steady state, with heating and cooling balancing each other. Fernandez 
& Reisenegger [6] calculated the simultaneous solution of f = i] npl = for a typical 
range of equations of state and found that, in a NS with a non- superfluid core, the pho- 
ton luminosity in the quasi-steady state depends only on the period and period derivative 
of the star, I* ~ (10 30 - 10 31 )(P_2o/j£s) 8/7 erg s _1 . Here, P_20 is the period derivative 
measured in units of 10~ 20 and P ms is the period in milliseconds. 
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FIGURE 1. Left panel: Thermal evolution with rotochemical (solid lines) and vortex-creep (long- 
dashed line) heating in classical pulsars. All curves correspond to stars with mass M = 1.4M©, magnetic 
field B = 2.5 x 10 n G, A18 + 8v + UIX* EOS, and initial temperature T = 10 n K. The abscissa 
correspond to the spin-down time (f s ), and the temperature upper limits are observational constraints for 
specific pulsars. Right panel: MSP with mass M = 1 .76M©, magnetic field 8 = 2.8x 10 8 G, initial period 
Pq = 1ms, and A18+5v+UTX* EOS. The solid and long-dashed curves correspond to rotochemical and 
vortex creep heating, respectively. The short-dashed curve shows the thermal evolution without heating 
(standard cooling). The error bar shows the temperature measured for the MSP J0437-4715. 

RESULTS AND DISCUSSION 

In order to solve the thermal balance equation (1) and generate the evolutionary curves, 
we use the computational code of Fernandez & Reisenegger [6]. This considers a 
realistic model of NS structure, with the core composed of neutrons, protons, electrons 
and muons. In order to model the rotational evolution, it is assumed that the rotational 
braking is due to the magnetic dipole radiation, without magnetic field decay. 

Thus, using an equation of state (EOS) that allows only modified Urea reactions, 
A18 + 8v + UIX* EOS, which makes the cooling relatively slow in comparison to 
EOSs that allow direct Urea reactions, Figure 1, left panel, shows the thermal evolution 
including vortex creep and rotochemical heating in classical pulsars. In the vortex creep 
mechanism, due to the uncertain value of the pinning energies for the interaction vortex- 
nuclei in the inner crust, the parameter J is adjusted to the temperature upper limit of 
pulsar B0950+08, wich constrains it to J < 2.8 x 10 43 erg s. As the luminosity of this 
mechanism depends only on the angular velocity derivative Cl, the resulting thermal 
evolution, for ages > 10 7 yr, is independent of the initial conditions and the previous 
thermal history. On other hand, in the case of rotochemical heating, the high magnetic 



field of classical pulsars, B > 10 G, causes strong braking of the pulsar, inducing high 
chemical imbalances in the beginning of the thermal evolution that later are slowly 
reduced by chemical reactions. Because of this, in the classical pulsar regime, the 
thermal evolution with rotochemical heating is very dependent on the uncertain initial 
rotation period Pq of NSs. 

Similarly, the right panel shows the thermal evolution for MSPs. As in the classical 
pulsars, the vortex creep heating is independent of initial conditions or previous thermal 
history. In this case, we assumed that the parameter J is universal. Thus, in the thermal 
evolution, we used the previous upper limit on J imposed by the pulsar B0950+08, 
resulting in a "phenomenological" temperature prediction very near that observed in the 
pulsar J0437-4715. On the other hand, for rotochemical heating, due to the relatively 
weak magnetic field, B ~ 10 8 G, the chemical imbalances induced by the rotational 
braking grow slowly, generating chemical reactions at high ages, t > 10 7 yr [6]. Through 
this, the star arrives at a quasi-steady state in which cooling and heating are balanced, 
and where the thermal evolution is independent of initial conditions and only depends 
on the current value of the product 

So, we find that the vortex creep and rotochemical heating can be important both for 
classical pulsars and MSPs. The temperature observed in the MSP J0437-4715 is slightly 
above the temperature prediction for both mechanisms. However, the temperature pre- 
dictions of rotochemical heating can be raised if a superfluid core is considered in the 
models [10]. Finally, better constraints on the temperature of some classical pulsars such 
as B0950+08 could rule out vortex creep heating as the generator of the thermal emission 
detected in the MSP J0437-47 15. 
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